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It has been shown that, for all dimensions and signatures, the most general first-order linear 
symmetry operators for the Dirac equation including interaction with Maxwell field in curved back- 
ground are given in terms of Killing- Yano (KY) forms. As a general gauge invariant condition it is 
found that among all KY-forms of the underlying (pseudo) Riemannian manifold, only those which 
Clifford commute with the Maxwell field take part in the symmetry operator. It is also proved that 
associated with each KY-form taking part in the symmetry operator, one can define a quadratic 
\ function of velocities which is a geodesic invariant as well as a constant of motion for the classi- 

cal trajectory. Some geometrical and physical implications of the existence of KY-forms are also 
C^) ' elucidated. 
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I. INTRODUCTION 



In many evolutions taking place in a flat or curved background, isometries of the underlying space-time metric 
lead to conservation laws that also have clear geometrical meanings expressed by means of their local generators, 
Killing vector fields. As space-time transformations, flows of these fields specify the conserved quantities as their 
flow invariants. However, since the beginning of the 1970s, it has been recognized that many interesting properties 
^Jy of a given space-time are intimately related to hidden symmetries of its metric, which make themselves manifest 
—t] in higher rank tensorial objects that also provide additional conservation laws [3, 0, 0]. The building blocks of the 
mathematical structure behind these hidden symmetries and associated conserved quantities are the symmetric Killing 
tensors, KY-forms, conformal Killing tensors and conformal KY (CKY)-forms. Their defining relations are natural 
generalizations of those of Killing vector fields and of conformal Killing vector fields. Although these higher rank 
| objects can be completely determined by the metric itself (see for instance [3, Q), they are not directly related to 
^C) . isometries or conformal transformations. But they comprise, by definition, their generators in the rank one sector 
of their hierarchy. In most of the cases these additional conserved quantities lead to complete integrability of the 
considered problem. This important step was initiated with the works of Penrose and his collaborators [H, They 
have shown that it is the existence of a second-rank symmetric Killing tensor which can be written as the "square" 
of a KY 2-form that leads to Carter's fourth constant of motion, which is responsible for the complete integrability of 
the geodesic problem in the Kerr geometry Jp . For other research fields utilizing KY and CKY-forms and for earlier 
| references for these forms we shall refer to pi || [t| E3] and the references therein. 

An important place where some or all of these tensorial objects enter the analysis is the study of symmetry operators 
for the Dirac-type equations describing the motion in curved background with or without additional interactions. By 
now it has become a well-established fact that while CKY-forms take part in symmetry operators, via the i?-commuting 
argument for the massless Dirac equation, KY-forms are indispensable in constructing first order symmetries of the 
massless as well as massive Dirac equation in a curved space-time. It is a four-dimensional Lorentzian space-time 
where most of applications have taken place. The first seminal studies in this context were carried out by Carter, 
McLenaghan and Kamran [111; El, EH • The results of these earlier studies, obtained in four dimensions for the massive 
or massless Dirac equation in the absence of electromagnetic interactions, were recently extended by Benn and his 
collaborators, to an arbitrary dimension and signature [3, EE Ell- 
in the absence of additional interactions, all KY or CKY-forms of the space-time take part in the symmetry operator 
without any extra restriction. However, when interactions are included, some additional conditions arise which restrict 
the possible forms that can enter into the symmetry operator. To the best of our knowledge, these restrictions for a 
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four dimensional curved background were found for the first time by McLenaghan and Spindel [17] . In searching for 
the most general symmetry operator commuting with the Dirac equation in the presence of an electromagnetic field, 
they found that the symmetry operator can be constructed from KY-forms of the underlying background, provided 
that they separately fulfill some conditions involving the field itself. These included some conditions found before 
by Carter and McLenaghan [l^l ■ Some of these conditions were also obtained before by Hughston et al in a slightly 
different context: in searching quadratic first integrals for the charged particle orbits in the charged Kerr background 



In this study we first show that the main results of McLenaghan and Spindel, that is the construction of the 
symmetry operators out of KY-forms, can be extended to an arbitrary dimension and signature. We then obtain a 
unified condition which allows one to specify which KY-forms can take part in the symmetry operator, and hence 
define hidden dynamical symmetry of the problem. This is an algebraic condition and can be stated as follows; a 
KY-form of the curved background enter the symmetry operator if and only if it Clifford commutes with the force 
(Maxwell) field. In particular, we solve all the consistency conditions and find a concise way to choose the gauge in 
order to make 0-form component of the symmetry operator constant. Owing to a non-integrated consistency condition, 
this point has remained ambiguous in the literature. 

Our results include the results of Benn and his collaborators when the electromagnetic field is turned off. Finally, 
we prove that the quadratic functions of velocities defined in terms of each KY p-form entering the symmetry operator 
is not only a geodesic invariant but also, for an arbitrary number of dimension and signature, a constant of motion 
for the classical trajectory. 

We shall mainly use the notation of [l8| and adopt the following conventions and terminology. The underlying 
base manifold is supposed to be an n— dimensional pseudo-Riemannian manifold with arbitrary signature. Covariant 
derivative of spinors, that is of sections of a bundle carrying an irreducible representation of the (real or complexified) 
Clifford algebra, with respect to the vector field X is denoted by Sx ■ Then the Dirac operator on spinors is S = e a Sx a , 
where the local co-frame {e a } is the dual to the tangent frame {V Q }. Summation convention over repeated indices 
will be used throughout the paper. Juxtaposing e a and Sx a , or any other operator, or form will denote the Clifford 
multiplication. When acting on forms Sx and 8 will be denoted, in terms of the pseudo-Riemannian connection V, 
as Vx and 



Here d = e°AVx a and S — —ix a ^x a denote the exterior derivative and the co-derivative written in terms of covariant 
derivative operator Vjf. A is the exterior product and ix will represent the interior derivative with respect to X 
whose action on an arbitrary p-form a is defined, for all vector fields Y), by 



For dual basis elements we have ix b e a = e a (Xb) = S£. 

The rest of the paper is organized as follows. In section II, the general form of the first order symmetry operator 
of the Dirac equation with a potential term is specified. This is achieved by constructing and then by solving all 
consistency equations except the equation for the 0-form component of the non-derivative term of the symmetry 
operator. The unified dynamical symmetry condition announced in the title is established in section III by analyzing 
higher degree components of the consistency equations. Special cases of this condition and the integration of the 
remaining 0-form component are also presented there. In section IV, the correspondence between KY and closed 
CKY-forms are studied and Yano vectors are introduced. Implications of the existence of Yano vectors related to 
symmetry analysis and to the global structure of the underlying space-time are given in the same section. In section 
V, the first integrals of geodesic equations, constants of motion of classical trajectories and their connection with the 
KY-forms and the mentioned dynamical symmetry condition are considered. Derivation of consistency equations and 
the contraction of curvature 2-forms with a Yano vector are given in two appendixes. In the Appendix B determination 
of upper bounds for the numbers of linearly independent KY-forms is also provided. Section VI concludes the paper. 
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(f-e°V Xa =d-S . 



(i x a)(Y 1 ,...,Y p _ 1 )=pa(X,Y 1 ,...,Y p _ 1 ) . 



II. FIRST ORDER SYMMETRY OPERATORS OF THE DIRAC EQUATION 



We set out to our analysis by considering the Dirac equation 



rmp 



(1) 



for a complex spinor field ip and propose the following first order linear symmetry operator 



L = 2u a S Xa + V , 



(2) 
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such that L Clifford commutes with $ + \A. Equation (1) describes the motion of a massive, charged and spin 1/2 
particle, with unit charge and mass m, interacting with the curved background encoded in the Dirac operator $ and 
force field represented by the potential form A. By fixing the charge from the outset, we assume the effective coupling 
to the spinor field. In the beginning A is allowed to be an arbitrary inhomogeneous form, and the consistency 
conditions are derived in this general context. Later on, A will be taken to be a 1-form. In that case, the term 
involving A in equation (1) describes coupling to the Maxwell field F = dA. A may also involve potential terms of 
non-electromagnetic origin such as those of conservative forces. In the latter case, F will be referred to as a force 
field. In equation (2) uj a and CI are Z2-homogeneous, both even or both odd forms which act, like A, on spinor fields 
by the Clifford multiplication. 

In recent studies 0, [H, EH it has been proved that use of the graded Clifford commutator [, ] considerably eases 
the calculations in analyzing the symmetries of the Dirac- type operators. For a p-form a and an inhomogeneous 
Clifford form (3 = Y] /3(g); this commutator is defined as 

K/?] = a/3-^(-l)%)«. 
<? 

If a is a one- form this transforms, in terms of the main involution rj (which leaves the even forms invariant and changes 
the sign of the odd forms) of the Clifford algebra, to [a, f3] — a(3 — f3 v a. Let us suppose that 

[&L]+i[A,L] = 0, (3) 

is satisfied and let us call L even (odd) when ui a and CI are both even (both odd). When L is even its graded 
commutator with the Dirac equation becomes the usual Clifford commutator [,]_ and irrespective of n, L is a 
symmetry operator. That is, it Clifford commutes with the Dirac equation and maps a solution to another. When 
L is odd it anti-commutes with Dirac equation and fails to be a symmetry operator. However, in such a case if n is 
even Lz is a symmetry operator since the volume form z anti-commutes with odd forms. (In such a case A must be 
an odd form.) 



A. The Main Consistency Equations 

To determine w a, s and CI in equation (2), we equate the symmetrized coefficients of the covariant derivatives of 
each order to zero in equation (3). The second order derivatives come only from the first bracket of (3) such that 

[$,L] = 2i Xb uj a {S 2 (X a ,X h ) + S 2 (X b ,X a )} + ■ ■ ■ , (4) 

where ix is the interior derivative, 

S 2 (X, Y) — SxSy — Sy x y , 

denotes the second covariant derivative and the ellipsis stands for lower order terms. Equating the coefficients of the 
second order derivatives to zero in equation (3) yields 

i Xb uJ a + ix°Ub = , (5) 

for all a, b = 1,2, ... ,n. This is satisfied if and only if uj a = ix°-w where tu is possibly a Z-inhomogeneous form. It has 
been shown in Appendix A that, in view of (5), by equating the coefficients of equal power of derivatives in equation 
(3) we obtain 

Vjfaw = i x *tp - i[A, ix-u] , (6) 
V 2 (A a ,A a V = {^ + 2i(i x .w)"V^ B i4-i[i4 ) n], (7) 

where <p = <jL> — CI. These constitute two main sets of the consistency conditions that specify the form of possible 
symmetry operators of the Dirac equation, for all dimensions and signatures, in which A is a general form. Henceforth, 
we take A to be a 1-form and write, for an arbitrary form a 

[A, a] — 2i^a . 

Here A is the metric dual of A, which in terms of the metric g of the background is defined by A(X) = g(A,X) for 
all X. 
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B. Solutions: Emergence of KY-Forms 

Let us first concentrate on equation (6), which for 1-form A can be rewritten as 

Vx«w = ix"{<fi + 2u A w) . (8) 

Applying ix a to both sides of this equation, we firstly see that Soj = 0, that is, u must be co-closed. On the other 
hand, by applying e a A we obtain 

duj = ir(tp + 2\i A uj) , (9) 

where the linear map it scales each form component by its degree : ir(a) = e a A ix a ct. The (p + l)-form component 
of equation (9) reads, for p = 0, 1, . . . , n — 2, as 

<P( P +i) = p~fY dUJ(p) ~ 2ii A w (p+^) ' ( 10 ) 

and for p = n — 1, as ip^ = n _1 cL;(„_i). In view of the last two equations, equation (8) implies that each p-form 
component of oj must obey 

v *^(p) = ^-j-ijcdwc.) , (11) 

which is the well-known KY-equation. 

From equation (10) and by the fact that uj is co-closed we also obtain 

p 

n {n) = (i - -)rfw ( „_i) , (13) 

n 

for p = 0, 1, . . . , n — 2. Note that for p = 0, 1, . . . , n — 1 we have 

*A n (p+i) = P V I w (p) • (14) 
These provide the non-derivative term, except for its 0-form component, of the symmetry operator in terms of KY- 
forms and the potential. 

We now turn to the implications of equation (7). Differentiating equation (8) once more, we obtain 

V 2 {X a ,X a )u = -S(ip + 2\i A uj) , (15) 
and by combining this with equation (7), we arrive at 

i A n - = s(i A uj) + {ix^TVxA . (16) 

This can be used to obtain the 0-form component of and possible relations among its higher degree components. 
For this purpose we apply the general relation 

[5,i x ]+ = -ix"iv Xa x , 

to KY-forms and obtain 

Wa^) = ~ix>iv Xa A u ■ ( 1? ) 

We now make use of 

(i X au) v K a = K a A i X aUJ - i~ K jx*u, (18) 

where K a 's are 1-forms. This relation is a direct result of the standard Clifford multiplication rule for a right multi- 
plication of an arbitrary form by a 1-form. Using (18), by taking n a = Vx a A, and (17) in (16) we get 

i A Sl- ^dfl = (V Xa A) Ai x «u) . (19) 

The 1-form component of this relation reads as 

dfi(o) = -2i(i^fi( 2 ) - ix«U(i)Vx a A) , 

= -2i(V^ (1) -V* (1) A). (20) 

In obtaining the second equality we made use of (14). The equation (20) will be integrated in the next section after 
the higher degree components of equation (19) are further analyzed to uncover the symmetry condition stated in the 
title. 
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III. A UNIFIED DYNAMICAL SYMMETRY CONDITION 

In this section we first show that equation (16) contains an important algebraic condition which plays a prominent 
role in deciding which KY-forms can take part in the symmetry operator. To see this, we first recall that £1 = dui — ip 
which implies 

dfl = —dip = 2\di A u> , 
i A Q = i A dio — i A ip — ni A ip . 

In the last equality we made use of (9). On substituting these relations into (16), we get 

n(i A tp) = -(fi A uj + (jx.w)"Vx a i • (21) 

We now make use of 

[<f,ix] + = Vx+e a i VxaX , (22) 

where [,] + denotes the Clifford anti-commutator. Using (18), once again by taking n a = Vx a A, and (22) in the 
equation (21) we arrive at 

n(i A ip) -i A dl) + V A w = ~[dA,u]- , (23) 

where we have also used the following relation: 

(K a M X a - e a Ai So )w = ~[dA,w]- . (24) 

It is now easy to verify that, by the KY-equation (11) and by the equation (10), the left hand side of (23) vanishes 
and we obtain an important condition [dA, = 0. On the other hand, as is apparent in equation (24), since the 
Clifford commutator of an arbitrary form by a 2-form does not change the degrees of its components, each p-form 
component of u must Clifford commute with the force field F — dA: 

[f> (p) ]_=0. (25) 

This gauge invariant condition means that among the KY-forms admitted by the underlying space-time, only those 
which satisfy the above condition take part in the symmetry operators. Moreover, it must be emphasized that it is 
the force-field, not the potential form itself, that plays a selective role in this regard. 

At this point we should note that, in view of equations (10), (14) and equation (21) derived above, the higher degree 
components of equation (19) lead us again to condition (25) and yield nothing new. 

In order to integrate equation (20) as well as to obtain more practical statements resulting from the equation (25), 
we rewrite it as 

ix«F Ai Xa LJ {p) = . (26) 

This is obviously satisfied for p — and also for p = n. Hence the extreme cases p — 0, n do not impose any condition 
on F and on possible KY-forms that take part in the symmetry operator L. In general, a KY p-form u>( p ) takes part 
in the symmetry operator if and only if F is in the kernel of the operator ix a ^>( p ) A ix a - More practical refinements 
of this condition are attained for intermediate values of p. We first should note that the 0-form component W( ) of 
a KY-form can be any function and ui^ is the dual of a Killing vector field. Moreover, W(„) is a constant (parallel), 
that is, it is a constant multiple of the volume form: W( n ) = kz. For p = 1 we can write u>^ = K, where if is a 
Killing vector field. In this case, equation (26) reduces to 

ikF = , (27) 

which implies that F remains invariant under the flow generated by K . That is, CkF = 0, where Ck denotes the Lie 
derivative in the direction of K. In other words, even the generator of the isometrics will take part in the symmetry 
operators if they fulfill condition (27). 

In terms of the potential 1-form A the condition (27) reads as %k (e a A V Xa A) = which, in view of i^ x kA = 
»lVxW(i), can be evaluated to obtain 



= \7 K A - di K A - V^W(i). 
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By comparing this relation with equation (20) and then by integrating we obtain 

Q = 2ii K A , (28) 

up to a constant. This relation, which was not recognized in the literature before, provides a concise way of choosing 
the gauge in order to make constant. In such a case, only duals of Killing vector fields whose flows preserve the 
potential can appear in the symmetry operator. An equivalent but more instructive way of deriving (28) goes as 
follows. Recalling the fact that for a 1-form j3 the metric dual of V x /3 is Vx/3, we can rewrite condition (20) as 

dh (0) = -2i[A,K] L = 2i£ K A . 

Here denotes the usual Lie bracket of vector fields, and we have also used the zero torsion condition. Since K 
is a Killing field, the above relation simply reads as d&ltm = 2iCkA and then the usual action of Lie derivative on 
differential forms produces the desired relation in view of (27). 

For a KY 2-form W( 2 ) = 2~ 1 uj a i 3 e a A e b the condition (26) reads as io a ^F^ = for all values of b and c where the 
square bracket stands for the anti-symmetrization of the enclosed indices and F c f, are the components of F in the 
same {e a } basis. This condition is the tensor-language version [12j of the condition first obtained by Hughston et al 
Our condition (25), or equivalently (26), are natural generalizations of this and similar conditions. This reflects 
the efficiency of the Clifford calculus in this context. 

To discuss condition (26) for higher forms it is convenient to rewrite it, in terms of (p + l)-forms defined by 

0a = F A ix a U( p ) , 

as ix a Pa = 0. This is also satisfied for p = 0, n. In the latter case f3 a — are identically satisfied, since all /3 Q 's are 
(n+ l)-forms. For p = n — 1 we have /3 a = f a z for some set of functions f a such that equation (26) amounts to iyz = 
for V = f a X a . As the map ip z defined by (p z (V) = iyz between the vector fields and (n— l)-forms is an isomorphism, 
iyz vanishes if and only if V = 0, that is if and only if all f a vanish. This is equivalent to F A ix a <-^( n -i) = for all 
X a , or to 

i Xs FAw (B _ 1) =0. (29) 

In the case of p — n — 2, /3 a 's are (n — l)-forms and each one can be written, in terms of uniquely determined 1-form 
Ca = 0a;,e b , as j3 a =* a a . In such a case equation (26) amounts to <j a A X a = 0, which implies that cr a {, = o\, a . This 
can succinctly be expressed as 

F A (X a A i Xb -X b A ix a )^( n -2) = , 

for a, b = 1, . . . , n. Note that the considered cases p = 0, 1, n — 2, n — 1, n exhaust all possible forms of equation (25) 
in four dimensions. More appealing versions of some of these conditions will appear in the next section (see equation 
(34) and the remarks followed). 

Finally in this section we should note that, when A is zero we recover exactly the results of Benn and Kress for all 
dimensions and signatures, obtained in [3, [H, [IB] . In such a case, since the condition (25) disappears, all KY-forms 
of the underlying space-time take part in the symmetry operators and all of the non-derivative terms explicitly given 
by equations (12) and (13), except Oq which is constant, are exact. 



IV. CORRESPONDENCE BETWEEN KY-FORMS AND CLOSED CKY-FORMS 



The appearance of higher rank KY-forms in the symmetry operator indicates the presence of dynamical symmetries 
which are not isometries. We shall now establish a general one-to-one correspondence between the KY-forms and 
closed CKY-forms and, as a particular case, we shall show that for each KY (n — l)-form there exists a uniquely 
determined conformal transformation. This transformation is generated by a locally gradient field whose integral 
curves are pre-geodesics. In this way we shall see that the higher KY-forms are related to special CKY-forms and, in 
the case of (n — l)-forms, to the special conformal transformations rather than isometry transformations. 

For the two points mentioned above, let us consider the CKY-equation 

VxP(p) = p~Tl lxdp ^ ~ n^p+l * A Sp (p) ' (30) 

which has the well-established Hodge duality invariance and conformal covariance. A p-form pi p \ is called CKY p-form 
if and only if it satisfies equation (30) for all vector fields X . It is obvious that a p-form is a KY p-form if and only if 
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it is a co-closed CKY p-form. A less obvious fact is that a (n — p)-form is a KY [n — p)-form if and only if it is the 
Hodge dual of a closed CKY p-form. Indeed, the Hodge dual of (30) for closed CKY forms yields 

Vx„ *P( P ) = r-f ix a d *P(p) ' 

v ' n — p + 1 v ; 

where we have made use of the relation ix *4> =* {4> A X) which can be considered as the definition of the Hodge 
map. Since p( p ) is closed it is locally exact such that = da, where the (p — l)-form a may be termed as the 
KY-potential for the KY (n — p)-form o>( n _ p ) =* da. Note that 0-forms of KY and CKY coincide and that any CKY 
n-form is of the form = fz, where / can be any differentiable function. As a result, the Hodge map establishes 
a vector space isomorphism between the vector space of all KY p-forms and that of all closed CKY (n — p)-forms for 
all values of p. 



A. KY (n — l)-forms and Yano Vectors 

In the case of KY (n — l)-forms we can write 

W( n _i) = iyz =* Y , (31) 

where the vector field Y is the metric dual of the associated closed CKY 1-form Y. Y will be referred to, following 
McLenaghan and Spindel, as the Yano vector. In that case, connection with the conformal transformations can be 
established in a more instructive manner. As the space of (n — l)-forms and that of 1-forms are isomorphic, exactly as 
in (31), any form of the former type can be written as the Hodge dual of a uniquely determined latter type. Recalling 
that the divergence of a vector field U with respect to the volume form z is defined by Ljjz — (div z ll)z, we can write 

Vx„W(„_ii = -(div z V)*X a 
n 

for the KY (n — l)-form u)r n -i) =* V (where V need not be a Yano vector). Equivalently, by the commutativity of 
the Hodge map and covariant derivative, we have 

Vx a V = -{div z V)X a . (32) 
n 

On the other hand, for an arbitrary vector field U, VxU can be decomposed as follows: 

VxU = hxdU--XSU + Tx(U) , (33) 
2 n 

(see [lH section 6.13 and 0]) where the 1-form Tx(U) is defined by 

r*(tf) = \[{L ug ) - -tr(L ug )g](X) . 
2 n 

Thus, for each vector field V corresponding to a KY (n — l)-form we have, in view of (32), ixdV — — 2Tx(V). That 
is, V is conformal if and only if V is closed. 

Let us return to equation (31): Y is closed, obeys the relation SY = —div z Y and generates conformal transformations 
Lyg = 2Ag with the conformal weight 

A = h tr{LY9) = ^ L y9){x\x a ) . 

Y is locally a gradient field for Y is closed. From (32), it also follows that 

V Y Y = -{div z Y)Y , 
n 

that is, the integral curves of Y are pre-geodesics and may be re-parameterized to become geodesies. If Y is divergence- 
free, then it is covariantly constant. We should note that in some literature, a curve whose velocity vector j satisfies 
Vyy = fy, for some function /, is termed geodesic, which we here call pre-geodesic by adhering to the nomenclature 
of [18|, For the Yano vector we have / = div z Y/n. 
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Substitution of (31) into (29) yields ix a F A* Y = 0, which in terms of components can be rewritten as 

F ab Y b = . (34) 

Note that condition (27) can also be rewritten as F ab K b — 0, where K b, s are the components of the Killing vector 
K. Thus, KY (n — l)-forms and 1-forms take part in the symmetry operator if and only if the corresponding Yano 
and Killing vectors are contained in the kernel of the matrix F ab . These imply that if F ab is non-singular, which 
can happen only in even dimensions, no KY (n — l)-form and 1-form can take part in the symmetry operator. Our 
condition for KY 2-forms can also be stated as: a KY 2-form will take part in the symmetry operator if and only if 
the corresponding anti-symmetric matrix (u) ab ) commutes with F ab . 



B. Contractions with a Yano Vector 



Contractions of curvature 2-forms R a b, Ricci 1-forms P a and conformal (Weyl) 2-forms C a b, defined for n > 3 by 
(below TZ denotes the scalar curvature) 

1 TZ 

Cab = Rab ~{P a A e b - P b A e a -e a A e b ) , 

n — 2 n — 1 

with a Yano vector enable one to reach decisive statements about the global structures of the underlying space-times. 
The mentioned contractions are found to be as follows: 

WRab = -[(iy-Pa)efc - (iyPb)e a } 

n — 1 

Y b i Y Pa = YaiyPb , (35) 

iyCab = t 7T7 7^[(R-g ac - P ac )eb 

(n — l)(n — 2) 

-(TZg bc - P bc )e a ]Y c + -^—(P a Y b - P b Y a ) , 
n — 2 

where P ac are the components of the Ricci tensor. The second relation easily follows from a second application of iy 
to the first relation and the last one also follows from the definition of C ab and the first relation of (35). A derivation 
of the first relation is given in Appendix B. Equations (35) are generalizations in the language of differential forms to 
an arbitrary number of dimensions and signature of the tensorial relations first found by McLenaghan and Spindel in 
the case of four dimensional Lorentzian space-times. 

If P ab = lZg ab /n such that TZ is constant, that is in Einstein spaces, we have iyC ab = 0. In four dimensions, the 
existence of a Yano vector on an Einstein space implies that the space is conformally flat or of Petrov type N. In higher 
dimensions, the equation iyC ab — for nonzero C ab is necessary but not sufficient condition for being Petrov type N. 
In such a case the space-time can also be of Petrov type II [22j . When Y is non-null we can write the second equation 
of (35) as iyP a = AY" a , where A = Y b iyP b /g(Y, Y). In such a case the Yano vector is an eigenvector of P ab with 
eigenvalue A = TZ/n. When Y is null, it is also null with respect to the Ricci tensor, in the sense that Y a iyP a = 0, 
in the directions of non-zero components of the Yano vector. In other words, in such directions contraction of Ricci 
forms with a null Yano vector projects it to an orthogonal direction. 



V. QUADRATIC GEODESIC INVARIANTS AND CONSTANTS OF MOTION 



The most important physical implication of the existence of KY-forms is the fact that they provide first integrals 
of the geodesic equation and may lead to quadratic functions of momenta that are invariant also for the classical 
trajectories. As is shown in the subsection B below, the condition (26) plays a crucial and unified role in establishment 
of the second point. 



A. First Integrals of Geodesic Equations 

The easiest way to see the first point mentioned above is to consider the relation [Vx, iy\ = ivxV which holds for 
two arbitrary vector fields X,Y. When X and Y are equal to a velocity field 7 of a geodesic, that is Vjj = 0, the 
covariant and interior derivatives commute and we obtain, in view of the defining relation (11), VjijLu = for any 
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KY (p + l)-form w. Therefore, the p-form a = i^ui and hence its "length" \a\ 2 defined by \a\ 2 — g p (a,a) remain 
constant along the geodesic 7. Here g p is the compatible metric induced by g in the space of p-forms and we assume 
go to simply multiply the two 0-forms (see [HI], Chapter 1). 

The constancy of \a\ 2 can be considered as a special case of a more general fact [20|. To see this, let us consider 
the symmetric bilinear form 

{X, Y) = g p (i x p, i Y p) = e* {i x P A* i Y 0) , (36) 

defined, in terms of a (p + l)-form j3 on the cartesian product of the space of vector fields. Here e denotes the sign of 
the determinant of g. In the case of a KY (p + l)-form u, one can easily verify the cyclic identity: 

V x A' (w) (y,Z) + SJ Y K {u>) {Z,X) + V z K {uj) {X,Y) = . 

That is, the symmetrized covariant derivatives of K^' vanish. This shows that to any KY (p + l)-form to is associated 
a symmetric bilinear form which is the Killing tensor generalizing the so-called Stackel-Killing tensor that 

corresponds to a KY 2-form, first recognized by Penrose and Floyd (the second tensor of (38) below). Since 

Vx [K {u>) {X, X)] = 2K {u>) (V X X, X) , (37) 

if("0 (<y ? 7) is constant along the geodesic 7. 

For KY 1-form oj^ = K a e a , 2-form = 2~ 1 uj ab e ab and n-form u>( n ) = where A: is a constant, the components 
of the corresponding Killing tensor are found, respectively, to be 

K a K b , -uj ac u c b , ek 2 g a b (38) 

In the case of KY {n — l)-form *Y, the components of Y ^ can be written in terms of the Yano vector Y as 

K ab Y) = 9n-2{ix a *Y,i Xb *Y)=eg 2 (Y AX a ,Y AX b ) , 

= e[g(Y,Y)g ab -Y a Y b ] . (39) 

This last Killing tensor, in particular, has remarkable properties worth mentioning. By definition Y ^(Y, Y) = 0, 
that is, Yano vectors are null with respect to the associated Killing tensor. When Y is non-null (that is, g(Y, Y) ^ 
), y ^ has a one dimensional kernel spanned by Y, and therefore it is singular. In that case every vector orthogonal 
to Y is an eigenvector with the eigenvalue eg(Y, Y) and the trace of Y ^ is e(n — l)g(Y, Y). Moreover, the properly 
normalized Killing tensor K' = Y ^ /cg(Y,Y) is an idempotent projector K' ab K' b c = K' ac having rank n — 1. When 

Y is null we have K^ b Y ^ = — eY a Y b and in such a case, it is a rank-1 nilpotent (K^^Kc Y ^ h = 0) projector, projecting 
vector fields to the direction of Yano vector. 



B. Constants of Motion for Classical Trajectories 

It is known that any rank-p symmetric (covariant) Killing tensor provides a degree-p polynomial of velocity that is 
a geodesic invariant [19(. But, the second rank symmetric Killing tensors that can be constructed from KY- forms, as 
in (36), have a distinguished property of providing a quadratic invariant of velocity along classical trajectories. In the 
remaining part of this section we shall prove this statement in the most general setting of this paper. To be precise, 
let us consider the quadratic function 

where u = C is the world velocity of a charged (charge assumed to be unit) material particle obeying the classical 
equation of motion 

V u u = —i^F . (41) 
m 

Here, u a — dx a /dr such that x a, s are the local coordinates of the world curve C parameterized by the proper time r 
and V u u represents the acceleration of the particle. Hence 

J:(/ (W) °C) = C*(d T )fM = V u [K^(u,u)] 

= 2K^{W u u,u) , (42) 
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and in view of (41) 

^(f^oQ = -u a u b F ac K^ c . (43) 
dr m 

What we are going to prove is the constancy of along the classical trajectory determined by (41). Before doing 
that, it would be illuminating to examine first some special cases. One can easily verify that the right hand side of 
(43) vanishes if the components of the Killing tensor given by (38) are used, provided that KY-forms employed in 

defining satisfy the condition (26). In the case of K ab Y ^ given by (39), the contraction with the Maxwell field is, 
in view of (34) 

F ac K^ )c = eg(Y,Y)F ab , (44) 

which is obviously anti-symmetric and also makes the right hand side of (43) vanish. So, as a particular result, all 
quadratic functions constructed from KY p-forms for p = 1, 2, n — 1, n are constant along the classical trajectory. 
These exhaust all possible cases in a four dimensional space-time and, to the best of our knowledge, these are all that 
can be found in the related literature. 

We shall now prove that is a constant of motion for the classical trajectories determined by (41) for all 
dimensions and signatures as well as for all KYp- forms uj^ obeying the symmetry condition (26). For the proof, we 
first take the Hodge dual of both sides of the equation (43) and write 

A(/M c)*l = -«°u 6 /„6, (45) 
dr m 

where the n-form I a b is defined as 

lab =* (Fa C K^ )c ) = ixJx a Fi X c(J A* ix b uj . (46) 

In passing to the second equality of equation (46) we have used (36). Thus, is constant along the world line C(r) 
if and only if the right hand side of the equation (45) vanishes. Evidently, the anti-symmetry condition I ao = —Iba is 
sufficient (and can also be shown to be necessary) for the right hand side of (45) to vanish. To show the anti-symmetry 
of I a b we rewrite it as follows: 

lab = -\ix a (ix c F M X cu) +ix c F Aixjxcu] A* i Xb u , 

= i Xc F A i x -(ix a uJ A* ix b u) - ix a u A i x F A i x - *ix b u; , (47) 
= —ix a u A ix c F A* (ix b w A e c ) . 

The first term in the square bracket of the first line vanishes because of condition (26), and the first term of the second 
line vanishes since it can be written as an interior derivative of a (n + l)-form. In the third line of (47) we have made 
use of the identity ix<= *ix b ^ =* (ix b w A e c ) which was also used in the previous section. In view of ix c F = F c ke k 
and of another Hodge identity a A* (i — (3 A* a which holds for any two p-forms a and /?, we obtain from (47) 

lab = -ix b uAF ck e c A* (i Xa ui A e k ) 
= i Xb LO A i Xk F A* (i Xa LU A e k ) . 

By comparing this with third line of (47) we obtain I a b = —I a b which was what to be demonstrated. 

An immediate corollary of the constancy of is, by virtue of (42), the relation K^{SI u u, u) = which can be 
interpreted as follows. The world velocity and acceleration of a charged material particle are perpendicular to each 
other, not only with respect to the metric, but also with respect to the symmetric Killing tensors associated with each 
KY-form satisfying the condition (26). 

VI. SUMMARY AND CONCLUSION 

In this study, the most general first-order linear symmetry operators of the Dirac equation including interaction 
with Maxwell field in curved background of arbitrary n-dimension and of signature are specified. We have shown 
that all coefficients forms w a 's of the symmetry operator L — 2ui a Sx a + ^ are given, in terms of an inhomogeneous 
KY-form ui, by uj a = ix°oj. The components of f2 are explicitly calculated and are given by the equations (12), (13) 
and (28). They depend on the exterior derivative of KY-forms and their contraction with the potential field A. We 
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have also found a unified, gauge invariant dynamical symmetry condition which states that among all the KY-forms 
of underlying curved background only those which Clifford commute with the Maxwell field can take part in L. When 
u} a and f2 are even, L itself, but when they are odd and n is even Lz is a first order symmetry operator which Clifford 
commutes with the Dirac equation and hence maps a solution to another. 

The special cases of the dynamical symmetry condition are also discussed as they may provide valuable insights in 
applications. In particular, the KY (n — l)-forms and 1-forms take part in the symmetry operator if and only if the 
corresponding Yano and Killing vectors belong to the kernel of the anti-symmetric matrix F a t corresponding to the 
components of the Maxwell field. These imply that if F a \, is non-singular, which can happen only in even dimensions, 
no KY (n — l)-forms and 1-forms can take part in the symmetry operator. Implications of the existence of Yano 
vectors in specifying the global structure of the curved background are also discussed. Finally it has been proved that, 
for all KY-forms obeying the dynamical symmetry condition, there exists a quadratic function of velocity (defined by 
the equation (40)) which is a constant of motion for the classical motion in an arbitrary dimension and signature. 

All of these results are expected to provide a unified framework for symmetry analysis and to serve as a firm base 
in studying symmetry algebras and related conserved quantities of the Dirac equation for a given specific curved 
background and force field. 
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APPENDIX A: CONSISTENCY CONDITIONS 

By direct calculations we obtain 

[# u] = 2Lu a S Xa +<JL} = L + tp , 
where Lu a — i X aLu. It proves convenient to take L as L — uS\ — ip, where Q = (foj — tp. Using the fact that 

the symmetry condition (3) transforms to 

[f,]- = [^l-i[A£]. (Al) 

TZ being the scalar curvature of the spinor connection, $ 2 acts on the spinor fields as (see [lH chapter 10) 

,f = S 2 (X a ,X a )-^K. (A2) 

Therefore on substituting the relations 

[$ 2 ,lu]_ = [S 2 (X a ,X a ),Lu]_ 

= 2V Xa u;S Xa +V 2 (X a ,X a )u J , 

[$,<p\ = 2i X aipS Xa +fo , (A3) 

[A,L] - 2[A,i x .u]S Xa -2(ix*u) v Vx.A+[A,Sl], 

into equation (Al) and then on equating the coefficients of equal power of derivatives, we obtain the consistency 
conditions (6) and (7) of the main text. 

APPENDIX B: CONTRACTION OF CURVATURES WITH A YANO VECTOR 

By differentiating the KY equation (11), the action of the Hessian (see the remark following equation (4)) on a KY 
p-form is found to be 

V 2 {X a ,X b )uj (p) = -^—i Xh V Xa duj [p) . (Bl) 
p + 1 

Since the difference of V 2 (X a , X&) and V 2 (Xi>, X a ) is the curvature operator R(X a , Xf,), from (Bl) we obtain 

R(X a ,X b )iu {p) = — —{i Xb V Xa -i Xa Vx h )du{py (B2) 

On the other hand, the action of the curvature operator on any form a is known to be (see [l8j equation (8.1.11) and 
0) 

R(X a , X b )a = -i X kR ab A i Xk a . (B3) 

Note that the action on the one form X simply is —i x R a b- Using (B3) in (B2) and then by multiplying the result 
with e a A, we obtain 

p + 1 

V Xa dio (p) = R b a A i Xb uj (p) . (B4) 



Using this in (B2) for W( n _i) =* Y we find 



R(X a ,X b )*Y = L^-[i Xb *Y AP a - i Xa *Y A P b ] (B5) 
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where we have used the contracted Bianchi identity %x h Pa = ix a Pb- Let us first note that 

R(X a ,X b )**Y = e(-l) n i Y R ab , (B6) 
*[*(YAe 6 )AP ] = i Pa **(YAe b ) , 

= €[{i Y Pa)e b -Yi Xb Pa] , (B7) 
where e is the sign of det g. In deriving the above relations, we have made use of (B3) and of the following identities: 

**a (p) = e(-l)^"-f)a (p) , i Xa *Y =* (Y A e a ) . (B8) 

If we now take the Hodge dual of both sides of (B5) and use (B6), (B7) we obtain the first relation of (35) of the main 
text: i Y Rab = (n - l) _1 [(i y P a )e fc - (iyP fc )e a ]. 

As an aside, it is worth mentioning that if (B4) is substituted into (Bl), we obtain 

V 2 (X a ,X b )uj {p) = U Xb {R C a^ixu {p) ), (B9) 

which shows that the second covariant derivative of any KY p-form is determined by curvature characteristics and 
the form itself. This implies that the value of a KY p-form at any point is entirely determined by the values of the 
form itself and its first covariant derivatives at the same point. These remarks can be used to determine the upper 
bound for the numbers of linearly independent KY p-forms to be the Binomial number C(n + l,p + 1) 0, l2lj |. 



